Photoemission spectra of LaMnOa controlled by orbital excitations 
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We investigate the spectral function of a liole moving in the orbital-ordered ferromagnetic planes 
of LaMnOa, and show that it depends critically on the type of orbital ordering. While the hole does 
not couple to the spin excitations, it interacts strongly with the excitations of Cg orbitals (orbitons), 
leading to new type of quasiparticles with a dispersion on the orbiton energy scale and with strongly 
enhanced mass and reduced weight. Therefore we predict a large redistribution of spectral weight 
with respect to the bands found in local density approximation (LDA) or in LDA+U. 



Transition metal (TM) compounds are well known for 
their large diversity and richness in phenomena . This 
variety of properties is not only due to the strongly cor- 
related nature of the electronic 3d states in these sys- 
tems, often rendering them magnetic, and to the strong 
hybridization with the extended ligand valence states, 
but also due to the orbital degeneracy of the open 3d 
shells. In a localized system such orbital degeneracy 
will be lifted in one way or another - this is the well 
known Jahn- Teller (JT) effect. In concentrated systems 
this often leads to structural phase transitions accom- 
panied by a certain ordering of occupied orbitals. Al- 
ternatively, strong correlations may lead to orbital order 
via superexchange interactions 1^, which may dominate 
the JT related contribution . Effects of this kind are 
observed in TM compounds with three-fold tg orbital de- 
generacy, for example in LiV02 Hj, but are strongest in 
compounds with a two-fold Cg orbital degeneracy, con- 
taining Mn3+ or GP+ (d^), Co^+ or (d^), or Cu'^+ 
(d^) ions in octahedral coordination. Some well-known 
examples are cubic (LaMnOa, KCrFs and KCuFa) and 
layered (LaSrMn207 and K2CUF4) perovskites ||,^. 

Superexchange interactions stabilize the magnetic and 
orbital order and the low-energy excitations of the 
system are magnons and orbital waves (orbitons). For 
simplicity we shall neglect the coupling to JT phonons. 
In fact, it has been argued by Dagotto et al. that if the 
Hund's rule interaction Jh is large, the JT coupling leads 
to the same orbital pattern as superexchange, and we as- 
sume that its effect can be absorbed in the renormalized 
parameters of the electronic model. 

In this Letter we present a study of single hole mo- 
tion in a ferromagnetic (FM) orbital-ordered plane. Due 
to strong correlations the hole motion is restricted by 
the magnetic order In the FM and orbital ordered 
state p the hole can move without disturbing the back- 
ground due to the off-diagonal hopping between eg or- 
bitals, and is essentially predicted by the band struc- 
ture theory using either the local density approximation 
(LDA), or its modification for the correlated oxides, the 
LDA+U method |lO|. The picture that arises from the 



present work is, however, quite different. In addition 
to free propagation, the hole motion may frustrate the 
orbital order, just as it does for the antiferromagnetic 
order in the t-J model. As we show below, the inter- 
action between the hole and the orbitals is so strong in 
the manganites that propagating holes are dressed with 
many orbital excitations and form orbital polarons that 
have large mass, small bandwidth and low quasiparticle 
(QP) weight 

Ground state — The ground state of LaMnOa is formed 
by FM layers which stagger along the c-direction, also 
known as A-type antiferromagnetic phase. As a con- 
sequence of the double-exchange the motion along 
the c-axis is hindered and propagation of holes occurs 
only within the FM (a, 6)-planes. Consider first the or- 
der of degenerate Cg orbitals induced by superexchange 
in a two-dimensional (2D) FM plane with one electron 
per site (n = 1). The orbital degrees of freedom are most 



conveniently described by pseudospins: 



and 



Tf — ^a-f , where erf are Pauli matrices, and we label 
the orbital \x) = \x'^ — y^) and \z) = |3z^ — r^) states 
by the pseudospin components: ||) = \x) and |J,) = \z). 
In the limit of large on-site Coulomb interaction U, the 
effective interaction between the orbital pseudospins is 
given by the superexchange J > 0, which involves the 
high-spin \^Ai) excited states |l3[. One finds 



(1) 



where (ij) denotes a nearest neighbor pair, and the nega- 
tive (positive) sign applies to a bond along the a(6)-axis. 
In the ground state the occupied eg orbitals are stag- 
gered: l^i) = i\x) + \z))/V^ and \iA = {\x) - \z))/V2 on 
A and B sublatice, respectively jl^. These orbitals de- 
fine locally a new basis obtained from {\x), \ z)} orbitals 
by a rotation by the same angle on the two sublattices, 

V'A = V-B = f 0. ^ 

We write the kinetic energy for holes introduced in the 
orbital-ordered state in terms of this new orbital basis, 
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FIG. 1. Electronic structure in the reduced BZ 

\X = (7r,0), S = (7r/2,7r/2)] as obtained from Ht (2) for: 
L'^ = (left), and U — 5t (right). In the bottom part 
the lower Hubbard band (centered at a; = 0) obtained for 
U = lot (dashed lines) is compared with the (7 — > oo disper- 
sion £^(0 = 0) (full lines). 



(2) 



where fj^ creates a hole on site i in orbital a, and up- 
per (lower) sign applies to a bond || a(6)-axis. Here the 
indices and 1 refer to the occupied (|0)) and unoccu- 
pied (|1)) orbital states, as determined by Hj. The hop- 
ping integrals in Eq. (^ follow from the Cg symmetry 
and therefore depend strongly on the pair of orbitals at 
nearest-neighbor sites involved in a hopping process. 

LDA + U bands — The Hamiltonian Ht gives a metallic 
system with a bandwidth of 6i, and can be seen as a tight- 
binding representation of the Cg bands in an idealized 
structure without lattice distortions. However, LaMnOs 
is a Mott-Hubbard system with large on-site Coulomb 
repulsion U which acts between the occupied and empty 
Cg states, Hu ^UJ2i fiofioflifn- I'^ the orbital-ordered 
state at n = 1 the charge fluctuations are suppressed and 
one may use the constraint //o/jq = 1. The Coulomb 
interaction acts then as a local potential on the unoccu- 
pied states, and the band structure is split into a lower 
and upper Hubbard band separated by a gap cx U, as 
is typical for a Mott-Hubbard insulator (Fig. |l|). For 
the manganites U/t « 10, and one finds that the lower 
Hubbard band has a dispersion ~ 2t, quite close to the 
U = oo limit. This resembles the strong redistribution 
of the bands in the LDA+U approach with respect to 
those found within the LDA, as obtained for LaMnOa 
[ p?! when the local electron-electron interaction terms 
cx U are treated in the mean-field approximation ]T8[ |. 

In our framework, however, the interpretation of Eq. 
(0) is quite different. On the site with the hole both 
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FIG. 2. Schematic representation of a hole in an or- 
bital-ordered ground state (a). The occupied (empty) orbitals 
1^) and Iv) are shown as filled (empty) rectangles. The hole 
can move either without disturbing the orbital order (b), or 
by creating orbital excitations (c). 

orbitals are empty [Fig. ^(a)]. An electron on a neigh- 
boring site can hop to either one of these orbitals. If the 
orbital order is preserved in the hopping process [Fig. 
0(b)], the hole propagates freely with a dispersion deter- 
mined by the first term /io/jo ^t- Unlike in the 
t-J model, the hopping couples the orthogonal states 
and on both sublatices. In contrast, if the orbital or- 
der is locally disturbed by the hopping process ~ fafjo 
[Fig. |^(c)], an orbital excitation occurs In this way 
the propagation of the hole is coupled to the orbital ex- 
citations determined by Hj. This coupling is a direct 
consequence of the strongly correlated nature of the eg 
electrons and its repercussions on the low energy electron 
dynamics cannot be taken into account within a mean- 
field (or LDA-I-U) treatment of the correlations. 
Orbital t-J model — Our total Hamiltonian, 



(3) 



includes, in addition to Hj (1) and Ht (2), a crystal-field 
term H^ = —E^ , which occurs, e.g., due to unaxial 

pressure. It acts like a magnetic field in the orbital sector 
which modifies the occupied orbitals and implies different 
basis rotations by V'a(b) = f T on A and B sublattice 
[ p^ , respectively, with sin 20 — —Ez/^J- 

Free hole dispersion — The hopping of the hole can be 
expressed by the change of the orbital background using 
the Schwinger boson representation [^ for the orbital 
variables which obey the constraint of no double occu- 
pancy in the limit oiU ^ oo: //q = ^io^I ^"^^ fli ~ ^ii^L 
with hi standing for a moving fermionic hole. In the or- 
bital ordered state the bio bosons are condensed, fo^o — Ij 
which leads to a free dispersion of the hole (Fig. |l|), 



(4) 



where e^('^) ~ t{—2sm2(f) + l)7k is determined by the 
orbital order, and = i(cos + cos ky). 

Orbital excitations — The orbital background is de- 
scribed by a local constraint for T = 1/2 pseudospins 
^io^io + ^Ii^ii — 2T, where 6Jq and are boson opera- 
tors which refer to the occupied and empty state at site i. 
The orbital excitations are calculated by using the lowest- 
order expansion of the constraint around the orbital- 
ordered ground state: ~ ^(^i + ^i) £^nd Tf ~ T — b\b^, 
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with bi = bii. The effective boson Hamiltonian is diago- 
nahzed by a Fourier and Bogohubov transformation de- 
fined by fcjj. — UyaiL+^Q^-k; with the coefficients {uk, 
given in Ref. Ilq]. One finds 



Ho = ^tJk(0)a[ak> 



(5) 



1 /2 

where LU\^{(p) = 3 J [l + i(2cos4(/)— l)7k] is the or- 
biton dispersion. This single mode, defined in the full 
Brillouin zone (BZ), is equivalent to two branches of or- 
bital excitations obtained in the folded zone in Ref. . 
The orbital excitations depend sensitively on the orbital 
splitting Ez- At orbital degeneracy {Ez = 0) one finds a 
maximum of Wk = 3 J(l 4- ^jkY^^ at the F = (0, 0) point 
and a weak dispersion ^ J. In contrast, for E^ — ±2J 
orbital excitations are dispersionless, with Wk — 3 J. 

Hole-orbiton coupling — The remaining part of Eq. (|^) 
describes the hole-orbiton interaction [Fig. §(c)] , 



ho 



k,q 



.^4.q"q 



^k.q^q+Q + ^-C] , (6) 



where Q = (7r,7r), A/k,q = 2 cos 2(/) (uq7k-q -l- Wq7k), 

^^k,q = -V3 (Uq77k-q - fq??k), and ?7fc,,fe„ = 7fex,fc«-Hir- So 

far we have reduced the total Hamiltonian (0) to an ef- 



fective Hamiltonian Ti., 



eff 



Hfi + Ho -\- Hfio, linearized in 



the slave fermion formalism, and treating orbital and hole 
dynamics on equal footing. Let us, for clarity, summa- 
rize the differences to the standard t-J model describing a 
hole in a quantum antiferromagnet: (i) the orbital model 
contains a free dispersion of the hole s'y.{(f>) which 
depends strongly on Ez; (ii) the orbital excitations are 
in general different when the momentum is changed by a 
nesting vector Q [a;k+Q(0) 7^ Wk(0)], and their dispersion 
varies with Ez , and finally, (iii) the hole scattering on the 
orbital excitations has a richer analytic structure than in 
the t-J model, with new processes cx iVk.q. A qualita- 
tively similar case to the t-J model arises at Ez = —2 J, 
where the free dispersion vanishes, but the orbital ex- 
citations are simultaneously dispersionless, so that the 
orbital model (||) reduces to the t-J^ model. 

We investigate the spectral function and QP properties 
using the self-consistent Born-approximation (SCBA) 
[ 21 1 , known to be very reliable for the single hole problem 
[20|. Treating TCeff in the SCBA, we find the selfenergy 

I](k, u;)^t^J2 [^'^k,k-q G(k - q, - cu^m 
q 

+ <k-qG(k-q,C.-C.q+Q(<^))] , (7) 

which, together with the Dyson equation for the hole 
Green function G(k, w), represents a closed set of equa- 
tions and was solved self-consistently by numerical iter- 
ation on a grid using 160 k-points. 

Quasiparticles — The spectral functions found in 
SCBA, A{'k,uj) = -ilmG(k,a;), consist of a QP band 



FIG. 3. Contour plot of the spectral function (left) and 
density of states DOS (right, blue lines) for J/t — 0.2, = 
(top) and Ez — 2 J (bottom). Quasiparticle states at the top 
of the valence band are highlighted in yellow color, while free 
hole bands £k(<?!') and DOS are shown by red lines. 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 

J/t 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 

J/t 



FIG. 4. Quasiparticle properties at l?^ = as functions 
of J/t: (a) inverse effective mass ra/m* (full circles) and the 
bandwidth W* /W (empty circles) in units of m and W de- 
fined by ek(0) (^; (b) weights of high-energy QP ar (filled 
circles) and low-energy QP bp (open squares) at the F point. 



close to the Fermi energy, while the excitations deep in 
the valence band are incoherent, taking typical values of 
J/t < 0.5. In Fig. ^ we present contour plots of the 
spectral function. The total width of the spectrum is 
comparable to the free dispersion, as shown by the den- 
sity of states, but the coupling to the orbitals leads to a 
strong redistribution of spectral weight. 

As the accurate orbital wave functions are not known 
we investigated the spectra for a few representative 
orbital ordered states. The QP part of the spectral func- 
tion obtained for the alternating (|a;) ± \z))/y/2 orbital 
order a,t Ez = resembles the free dispersion with a 
maximum at the F point, but its bandwidth is reduced 
to ~ J. When Ez = 2 J, which corresponds to the alter- 
nation of 3a;^ — r^/3?/^ — orbitals ||2^, the spectral func- 
tion changes markedly (Fig. |^). The incoherent weight 
is now distributed more smoothly and the original free 
dispersion Skl'?^) {^'^'i lines) can still be recognized in the 
spectral function but is strongly damped. Most strik- 
ingly, there is only one QP band with appreciable weight 
left in a limited region of the BZ, while the second band 
is absorbed by the continuum. In the opposite case of 



-2J, i.e., a;2 - z^/y'^ 



orbitals alternate, the 



spectrum (not shown) is identical to the ladder spectrum 
of the t-J^ model; it is k-independent and consists of a 
set of (5-functions at approximately equal energy intervals 
In this respect the spectrum found at Ez — can 
be viewed as a compromise between the extremes of the 
ladder spectrum, from which it retains some character 
of the enhanced/reduced spectral weight at regular en- 
ergy intervals, and the Ez = 2J spectrum in which the 
dispersive features are smeared out over the entire band. 

Finally, we analyze the QP properties at orbital de- 
generacy (Ez = 0). Two QP states at the F point 
have the spectral weights ar and br, respectively, and 
determine the bandwidth W* of the QP band (Fig. 
H). At low doping the states at the top of the QP 
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band are filled by holes, and thus the transport prop- 
erties depend on the effective mass found at the F point, 
(9^i?k/9A:^)|k=o oc m/m*, where i?k is the QP energy. 
In the weak-coupling regime [J > t/2) the QP bandmass 
is almost unrenormalized, while it is strongly enhanced 
in the strong-coupling regime (J < t/2). For the man- 
ganites J/t ~ 0.1 @, so that W* /W ~ 0.1, and the QP 
mass TO* is increased by an order of magnitude due to the 
dressing of a hole by the orbital excitations. The weights 
of the two QP branches at the F point are strongly re- 
duced in the strong-coupling regime, but the QP at the 
top of the valence band is still quite distinct. 

Summary — In conclusion, we have shown that a hole 
moving in FM planes of LaMnOs couples strongly to or- 
bital excitations which results in a large redistribution of 
spectral weight compared to single-electron or mean-field 
treatments. The coherent quasiparticle band and the in- 
coherent part of the spectral function depend critically 
on the orbital ordering and would therefore be different 
for FM planes of cubic and layered manganites. This 
prediction could be verified by angle resolved photoe- 
mission. Unfortunately, such experiments have not yet 
been performed in the orbital-ordered phase at low dop- 
ing concentration. Recent experiments by Dessau et al. 
psf on the highly doped layered system Lai.2Sri.8Mn207 
showed strong incoherent features which were discussed 
in terms of small lattice polarons. Our study points out 
that the incoherence observed in these experiments might 
be attributed instead to orbital excitations. 
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